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Abstract 

In the mean field approximation, we evaluate the temperature dependence 
of the anchoring energy strength of a nematic liquid crystal in contact with a 
solid substrate due to thermal fluctuations. Our study is limited to the weak 
anchoring case, where the microscopic surface energy is small with respect 
to the mean field energy due to the nematic phase. We assume furthermore 
that the physical properties of the substrate can be considered temperature 
independent in the range of the nematic phase. According to the thermo- 
dynamical perturbative approach, the macroscopic surface energy is deduced 
by averaging the microscopic one, with a density matrix containing only the 
nematic mean field. We show that the thermal renormalization of the an- 
choring energy coefficients is proportional to the generalized nematic order 
parameters. Our analysis shows also that the thermal renormalization of the 
anchoring energy coefficients predicted by means of Landau-like theories is a 
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first and rather rough approximation in the whole nematic temperature range. 
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One of the least-understood areas of physics and chemistry of hquid crystals concern 
the anchoring phenomenon and the temperature surface transitions at the interface liquid 
crystal-solid or soft substrate. These phenomena are important also from a practical point 
of view, since they play fundamental role in the realization of displays. In this Letter we 
analyse these phenomena in nematic liquid crystal media. 

Nematic liquid crystals are anisotropic fluids made by anisometric molecules having 
quadrupolar symmetry. Their intermolecular interaction, Vn, is such to orient the molecular 
axes, u, along a common direction n, called director . It coincides with the optical axis of 
the medium. When a nematic liquid crystal is in contact with a substrate, the orientation 
of n at the surface results from a balance of the anisotropic interactions with the bulk and 
with the substrate. In the absence of bulk distortions, the surface orientation of the nematic 
director coincides with the "easy axis", no 0. It is such to minimize the anisotropic part 
of the surface energy characterizing the interface between the nematic and the substrate. 
Long ago Bouchiat and Langevin-Cruchon found a strong temperature dependence of the 
easy axis. The measurements of Ref. |^ have been repeated by other groups with similar 
results Several models have been proposed to interpret this phenomenon. According 

to Parson [^] the easy axis results from the competition between dipolar and quadrupolar 
interactions, which depend on the temperature in different manner. In special situations a 
surface nematic orientation temperature dependent can be observed. The idea of Parson has 
been generalized by Sluckin and Poniewierski |TT|-[T3| and by Sen and Sullivan [|I4|. In all 



the models, the temperature surface transitions are due to a temperature dependence of the 
anisotropic part of the anchoring energy, which depends on the symmetry of the substrate 
and on the symmetry of the nematic phase. 

The aim of our paper is to analyze the temperature dependence of the anisotropic part 
of the surface energy. We assume that the nematic liquid crystal is not polar. From the 
molecular point of view it has a quadrupolar symmetry, whose principal axis coincides with 
u. The elements of the relevant tensor are qij = {3/2)[uiUj — (l/3)5jj]. In our analysis, we 
expand the surface energy in series of spherical harmonic functions. The coefficients of the 



expansion are the experimentally detectable anchoring coefficients. According to our model, 
all the anchoring coefficients of the same order depend on the temperature in the same 
manner. From this result it follows that in nematic liquid crystals the alignment transitions 
driven by the surface (the so called temperature surface transitions) are due to a surface 
anchoring energy which contains contributions of different orders. 

We analyze the temperature dependence of the anchoring energy using an approach based 
on the mean field theory. In our analysis we neglect all the inhomogeneities. We assume, 
furthermore that the surface potential is short range. 

Let us consider a surface molecule of the nematic liquid crystal. It is submitted to the 
mean field due to the other nematic molecules, whose corresponding energy is Vn, and to 
the interaction with the substrate, Vg. In this framework, the total energy, V, of a given 
surface molecule is = Vat + V5. If ~ Vs the extrapolation length b = K/W ~ aV^/Vs, 
where K is an average elastic constant and a a molecular dimension, is of the order of a 
molecular dimension [|T|. In this case, in the continuum limit it is possible to put 6 = 0, 
and assume that the surface nematic orientation is fixed by the surface interaction. This 
situation is known as the strong anchoring case, and it is not interesting for us here. The 
interesting case is the one in which Vn ^ Vs, corresponding to a situation where b ^ a. This 
case corresponds to the weak anchoring situation, to which we will limit our investigation. 
In our analysis the small parameter used to expand the surface energy in power series is 
Vs/^N ^ 1, in the weak anchoring situation. On the contrary, the surface scalar order 
parameter S is not supposed to be a small quantity. 

Vn describes the tendency of u, which defines the molecular orientation, to be oriented 
along the nematic director n. Usually, it is approximated by means of the Maier-Saupe's 
mean field |T^, V^^, according to which V^ oc riiqijUj = P2(n ■ u), where P2 is the second 
order Legendre Polynomial. In this framework V^ = —vP2{'n ■ u)S, where f is a molecular 
constant and 5* = (-P2(n • u)) the nematic scalar parameter. A generalization of the Maier- 
Saupe theory has been proposed by Humphries et al. According to this generalized 

mean field theory the nematic mean field is given by 
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VAf(n ■ u) 



(1) 



where V21 are molecular parameters, and 5*2; = {P2i{i^ ■ u)) the nematic order parameters, 
given by the self-consistent equations 



The Maier-Saupe potential, , is obtained from Vat putting V2k = v6i^k- 

The interaction connected to Vs has to describe the tendency of the surface to orient 
the surface nematic molecules along the "easy direction", Hq. This direction depends on the 
symmetry of the surface and on the molecular properties of the mesophase. Since we limit our 
analysis to non polar media, Vs has to be an even function of u. It follows that Vs is, actually, 
a function of the tensor 'q and can be written, in general, as Vs{u) = Vsiji) = J2k 'Wk{0)Li:{*q), 
where L^i^q) indicate the scalar quantities we can build with the molecular tensor of elements 
Qij = {3/2)[uiUj — {l/3)6ij] and the elements of symmetry characterizing the surface. Each 
term of the expansion of Vsiji) represents a given interaction, like induced dipole-induced 
dipole or quadrupole-quadrupole and so on ||18|; the "intrinsic" anchoring coefficients Wk{0) 
are physical parameters connected with the type of interaction described by Lk(q). Wk{0) 
refer to specific fundamental interactions, and are assumed to be temperature independent. 
In this case thermal effects arise only from the temperature dependence of the degree of 
alignment of the nematic molecules. This conclusion is valid only if in the temperature range 
of the nematic phase the physical properties of the substrate can be considered constant. In 
the opposite case Wk{0) depend also on the temperature, via the substrate. Since we assume 
that Wk{0) are temperature independent, our theory works well when the substrate is a solid 
crystal. Deviations from our prediction are expected for nematic samples oriented by means 
of surfactants. 

For our future considerations it is useful to describe the molecular direction and the 
nematic director in terms of the polar angles with respect to a cartesian reference frame 
having the z-axis parallel to the geometrical normal to the flat surface and the x-axis along 



S21 



Jq P2i{n ■ u) exp[-/3Viv(n ■ u)]d{n ■ u) 
Jq exp[-/5VAr(n ■ u)]rf(n • u) 



(2) 
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the possible surface anisotropy. Let 0, $ and 6*, be the polar and azimuthal angles defining 
u and n, respectively. Consequently 

Vs{q) = Vs{Q,<^) =Y.Wk{0)Lk{Q,<l>). (3) 

k 

By decomposing the functions 1/^(0, $) in series of spherical harmonics functions Y^{Q, $) 
we obtain Lfc(0, $) = X^m O/T^fc^l©) '^)- Since Lk = Lk(q) and hence Lfc(0,$) = Lk{Tc - 
0,7? + $) for all k, we deduce that k = 21. It follows that for non-polar nematic liquid 
crystals ^2^(0, $) = Hm '^ii^'M i.^ ) ^'^^ the microscopic surface energy can be written as 

^5(0,$) = E^2Ko)E«^>^2?(0,*)- (4) 

I rn 

The macroscopic anchoring energy 14^ (n) = W{9, 0) is obtained by averaging Vs over the 
molecular orientations u, or over and $. Since in the problem under consideration Vs ^ 
Vat, Vs can be treated as a perturbation. According to the thermodynamic perturbation 
theory [|T9| we have W{0, 0) = (^5(0, $)), and hence, as it follows from Eq.(|^), 

w{e, 0) = E ^2/(0) E <i{y^{^. (5) 

I m 

where {A) = Tr{pA)/Tr{p), and p = exp(— /JVat) is the density matrix. In order to derive 
the macroscopic surface energy W{6, 0) we have first to express Vs{Q, $) in terms of a polar 
coordinates system based on the director n as polar axis. The cartesian reference frame has 
to be rotated in such a way that z ' = n. We will indicate with -d, (f the polar angles of u 
with respect to the rotated coordinate system. In this case 

Yne,^) = j:DL,m'iO,<P)Yr\^,^), (6) 

m' 

where D^^/(6',0) are the elements of Wigner's matrix. Since (Y[^' {•& , ip)) = {Y^^{'d))5m'fl 
we obtain from Eq.(D (Y,™(0,<l>)) = o(^, 0)(Y;°(?9)). By taking into account that ^ 
L)^q(^,0) = r;"(^,0), we have finally, as it follows from Eqs.(||,|J|), 

w{e,(i)) = Y.W2i{o)SMe,<p), (7) 
I 



where we have taken into account that Y2i{'d,ip) = P2i{cos'd). Eq.(0) is a consequence of 
the fact that we regard all anisotropic effects as perturbation, so that they do not need to 
be included in the computation of the averages values. There is axial symmetry about the 
direction of n in the imperturbed system and only the member m = of the is different 
from zero. By comparing Eq.(|^) with Eq.(|^) we deduce that the temperature dependence 
of the parameters describing the anisotropic part of the surface energy is given by 

W2liT) = W2liO)S2l. (8) 

This means that the temperature dependence of W2i{T) / W2i{0) coincides with the tempera- 
ture dependence of the 2/-th scalar order parameter. 

According to the analysis presented above, where the macroscopic anchoring energy is 
given by the series expansion in spherical harmonic functions shown in Eq.(^, the thermal 
renormalization of the anchoring coefficients is given by Eq.(P). /^From these results it follows 
that the anchoring coefficients of the same order in the expansion have the same tempera- 
ture dependence. Consequently, in the frame of our model, temperature surface transitions 
are possible only in nematic samples whose anchoring energy contains contributions from 
different order in the spherical harmonic functions expansion. 

The ratios S21/ S vs. S, for / = 2, 3 and 4, in the Maier-Saupe approximation, can be 
easily evaluated in the nematic phase, where 0.4 < S < 0.8. A direct calculation shows that 
S21/S < 0.2, for / = 3,4, as it is shown in Fig.l. This explains why, usually, the anisotropic 



part of the surface anchoring energy given by Eq. (0) is well approximated by few terms ^1 

In the low temperature region, where —(3Vn ^ 1, the fluctuations of u with respect 
to n are small. In this region n ■ u = cos^ ^ 1 _ (l/2)'i9^ + (9(4), i.e. "i? <^ 1, and 
P2/(cos'/9) = 1 — [l{2l + l)/2]'i9^ + C(4). Consequently, from Eq.(^, the order parameter 5*2^ 
is found to be 

/(2/ + 1) f /(2/ + 1)1 

S21-I- ~ j-^-^l ' 

where B = (3J2k k{2k + l)v2kS2k- The main nematic scalar order parameter S = (-P2(n ■ u)) 
is given by S* = exp(3/i?),as it follows from Eq.(P). The other order parameters can be 



determined in terms of 5* by 5*2^ = 5''*^^'+^^/^. In the low temperature region, the thermal 
renormalization of the anchoring coefficient is then given by 

^2^(2^) =^^2^(0)^^. (10) 

In particular, in this range of temperature, W2{T) / W2{Q) = S and W4{T)/w/i{0) = S^^^^. 
The temperature dependence given by Eq.(|l^) reminds the Akulov-Zener law for magnetic 
anisotropy, well known in ferromagnetism theory p2| , pB| 

As an example, we consider now a nematic liquid crystal limited by an isotropic substrate. 
In this case only the polar angle 6 enters in the description. The analysis of the temperature 
surface transitions in a system of this kind is usually performed by means of a Landau's 
expansion of the anisotropic part of the surface energy |]TI| , ^ . According to this approach 
W^(n) is expanded in power series of the invariants made with the elements of symmetry 
characterizing the nematic phase (which is the nematic tensor order parameter of elements 
Qij = {3/2)S[ninj — (l/3)(5jj]), and the substrate (which is the geometrical normal z). 
In the Landau-like approaches the quantity playing the role of expansion parameter is S. 
However, since the nematic-isotropic phase transition is first order, S is never very small 
(at the transition point it is of the order of 0.3 At the second order in S, W^{6) = 



W2-P2(cos6') + w^P4{cos6) + 0{3), where = ^iS + 026*^ and = a^S'^, in which oi, 02 
and 03 are constant parameters, temperature independent |TT]JT^,|2^ . 



Now we want to compare the prediction of a Landau's expansion up to the second order 
in S with the result of our mean field analysis. In the case under consideration the angular 
functions ^2/(6', 0) reduce to L2i{0) = P2i{cos9), and 

w{e) = J2w2i{T)P2i{cose). (11) 
I 

/^From Eq.(^ we obtain W2{T)/w2{0) = S. This means that at the first order in S the tem- 
perature dependence of the anchoring energy deduced by means of symmetry considerations, 
W2, and by means of the mean field agree. However, for / = 2 there is a discrepancy between 
the two approaches. In fact, according to the mean field we have W2(T)/w2(0) = S, and 
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W4,(T)/w4{0) = 5*4 7^ 5*^, whereas the Landau's approach predicts the temperature depen- 
dencies W2(T) = aiS+a2S^ and w^(T) = a-^S"^. More precisely, it predicts a renormahzation 
of the coefficient of P2{cosd), by means of a contribution, and a temperature dependence 
of the coefficient of P4(cos^) hke S"^. Of course, in the hmit of small 5* the two predictions 
agree. In fact, if S* ^ 1 the renormahzation of P2{cos9) in can be neglected with respect 
to the linear term in S. Furthermore, in this approximation, S4 oc S"^. However, in the case 
of large S the discrepancy between the two approaches can be large. In the low temper- 
ature region, where it is possible to use the approximate expressions given by Eq. ([lO|) for 
the thermal renormahzation of the anchoring coefficients, our mean field approach predicts 
W2(T) oc S and w^iT) oc 5^°/^. In Fig. 2 we show S"^, predicted by Landau-like models, and 
^10/3^ predicted by our mean field theory in the low temperature region, vs. S. As it is evi- 
dent from this figure, our theory represents an improvement with respect to the Landau-like 
approaches in the whole temperature range. 

To conclude we stress the main results reported in the paper. We have shown that the 
renormahzation due to the thermal fiuctuations of the anchoring coefficients W21 is of the 
kind W2i{T)/w2i{0) = S21 where S21 is the 2/-th scalar order parameter. In the particular case 
in which the nematic phase is described by the Maier-Saupe theory, W2i(T) / W2i{0) coincides 
with the average value of the 2/-th Legendre polynomial. We have also shown that only 
at the lowest order in the scalar order parameter the simple approach based just on the 
symmetry of the problem agrees with our mean field approach. This is a consequence of 
the hypothesis of small S, over which is based the validity of the Landau-like expansions 
of iy(n) in power of S. We have proposed also approximate expressions for the thermal 
renormahzation of the anchoring coefficients, valid in the low temperature region, where the 
fiuctuations of the molecular directions with respect to the nematic director are small. 
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FIGURES 

FIG. 1. S21/S vs.S according to Maier-Saupe theory. In the usual nematic range, where 
0.4 < 5 < 0.8, Sq/ S, Ss/ S < 0.2. This explains why, usually, two terms are enough to ap- 
proximate the macroscopic surface energy. The contribution from higher harmonics disappears as 
a result of the thermal fluctuations of the nematic molecules with respect to the director. 

FIG. 2. and 5^°/^ vs.S. According to Landau-like models at the second order in S, the 
thermal renormalization of the surface energy in P4{cos9) is proportional to S'^. According to our 
mean field theory, it is proportional to S4. In the low temperature region, where the fluctuations 
of the molecular orientation with respect to the nematic director are smaU, S4 ~ 5^°/^ The 
figure shows that our approximate expression for 5*4 represents an improvement with respect to 
the Landau-like analyses. 
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